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Abstract 

We extend a previous analysis of flavor-changing fermion-graviton vertices, by adding 
the one-loop SM corrections to the flavor diagonal fermion-graviton interactions. Explicit 
analytical expressions taking into account fermion masses for the on-shell form factors are 
computed and presented. The infrared safety of the fermion-graviton vertices against ra- 
diative corrections of soft photons and gluons is proved, by extending the ordinary infrared 
cancellation mechanism between real and virtual emissions to the gravity case. These re- 
sults can be easily generalized to fermion couplings with massive gravitons, graviscalar, 
and dilaton fields, with potential phenomenological implications to new physics scenarios 
with low gravity scale. 
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1 Introduction 



The investigation of the perturbative couplings of ordinary field theories to gravitational back- 
grounds has received a certain attention in the past [U El EJ 0]. While the smallness of the 
gravitational coupling may shed some doubts on the practical relevance of such corrections, with 
the advent of models on large (universal [SEllT] and warped [3 [9]) extra dimensions and, more 
in general, of models with a low gravity scale JTU] their case has found a new and widespread 
support. This renewed interest covers both theoretical and phenomenological aspects that have 
not been investigated in the past. They could play, for instance, a significant role in addressing 
issues such as the universality of the gravitational coupling to matter [11], in connection with 
the Lagrangian of the Standard Model. On the more formal side, as pointed out in some stud- 
ies [121 [IHl HH [El [16] , the structure of the effective action, accounting for anomaly mediation 
between the Standard Model and gravity, shows, in its perturbative expansion, the appearance 
of new effective scalar degrees of freedom of dilaton type. These aspects went unnoticed before 
and, for instance, could be significant in a cosmological context. For this reason, we believe, 
they require further consideration. 

All these perturbative analyzes are usually performed at the leading order in the gravita- 
tional coupling. This is due to the rather involved expression of the operator responsible for 
such a coupling at classical and hence at quantum level: the energy-momentum tensor (EMT) 
of the Standard Model. Its expression in the electroweak (EW) theory is, indeed, very lengthy, 
and the classification of the several amplitudes (and form factors) in which it appears - at 
leading order in the EW expansion - requires a considerable effort. 

The previous vertices discussed in the literature have been the graviton-fermion-antifermion 
vertex (Tff) and the graviton-gauge-gauge vertex (TJJ'), with J being a generic neutral 
current. A general discussion of all the amplitudes and the renormalization properties of these 
vertices has been presented in the TJJ' case in [17J, while the Tff case has been analyzed, 
for fermions of different flavors, in [11]. This second work also contains a detailed discussion of 
the direct implications of the result for massive gravity, leading, in this respect, to interesting 
noteworthy conclusions. In particular, it was shown that flavor-changing interactions coupled 
to a gravity background are local only if the graviton is strictly massless, with a long-range 
contribution appearing in the Newton potential only if a graviton has a small mass. 

In the current work we are going to revive this program and extend this previous analysis 
of the Tff vertex pj] to the flavor diagonal case, presenting the explicit results for the EW 
and strong corrections to the corresponding form factors. These are given in a basis which 
partly overlaps with the previous tensor basis of the flavor-changing case [11] , but with some 
new additions. 
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We anticipate that these results are relevant for some phenomenological consequences that 
may affect studies on these interactions. One of the special contexts, which we plan to address 
in the future, is represented by the neutrino sector. Indeed, the analysis of the effects of the 
one-loop EW radiative corrections to the gravitational interactions of neutrinos is missing in 
the literature. Although there are previous studies dealing with more general structures for 
the neutrino-graviton interactions beyond the tree-level approximation, see [18j and references 
therein, none of these works takes into account the exact neutrino-graviton vertex at one-loop 
level in the SM, including the neutrino mass dependence. In particular, these corrections are 
expected to play some role in the neutrino physics in the presence of a gravitational background, 
mainly due to the new parity- violating structures induced by EW corrections. In this respect, 
we are planning to analyze in a forthcoming paper, the impact of the one-loop EW radiative 
corrections on the gravitational interactions of neutrinos, including the effects of the EW con- 
tributions to the off-diagonal flavor transitions [HI [19] . The inclusion of these corrections in 
the fermion-graviton vertex might also help in clarifying whether Majorana or Dirac fermions 
behave differently with respect to the gravitational interactions, a much debated topic in the 
litarature [18]. Moreover, it is also left open the possibility of extending our analysis to more 
general gravitational backgrounds, with the inclusion of a dilaton field. 

Our work is organized as follows. We briefly discuss the structure of the embedding of 
the Standard Model Lagrangian in a curved space-time, assuming as a background metric the 
usual 4-dimensional one. The discussion is rather general, and remains valid also for more 
general cases, which include dilaton backgrounds obtained from the compactification of higher 
dimensional metrics. Then we turn in section [3] to illustrate the structure of the perturbative 
expansion, organized in terms of the various contributions to the EMT of the Standard Model. 
These are separated with respect to the particles running in the loops, which are the photon, 
the W and Z gauge bosons and the Higgs field. We conclude this section with a classification 
of the relevant form factors. 

In section H] we briefly discuss the derivation of an important Ward identity involving the 
effective action which is crucial to test the correctness of our results and to secure their consis- 
tency. Section [5] addresses the issue of the renormalization of the theory, which complements 
the analysis of |20j . We recall that no new counterterms are needed - except for those of the 
Standard Model Lagrangian - to carry the perturbative expansion of EMT insertions on cor- 
relators of the Standard Model, under a certain condition. This condition requires that the 
non-minimal coupling (x) of the Higgs field be fixed at the conformal value 1/6. We then 
proceed in section |6] with a description of the expression of the form factors for each separate 
gauge/Higgs contribution in the loop corrections. 

In section [7] we give a simple proof - at leading order in the gauge couplings - of the 
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infrared finiteness of these loop corrections, once they are combined with the corresponding real 
emissions of massless gauge bosons, integrated over phase space. The proof is a generalization 
of the ordinary cancellation between real and virtual emissions, in inclusive cross sections, to 
the gravity case. Finally, in section [8] we give our conclusions. 



2 Theoretical framework 

We recall that the dynamics of the Standard Model plus gravity is described by the Lagrangian 

S = S G + S SM + Sj = ~ J d 4 x^R + J d A x^g~C S M + xj d 4 x^RH^H , (1) 

This includes, beside the Einstein term So where R is the Ricci scalar, the Ssm action and a 
Si term involving the Higgs doublet H. The latter is responsible for generating a symmetric 
and traceless energy-momentum tensor [21], also called "term of improvement". Ssm, instead, 
is obtained by extending the ordinary Lagrangian of the Standard Model to a curved metric 
background. 

Notice that Sj vanishes in the flat space time limit, due to the vanishing of the Ricci scalar 
in the same limit, but the EMT which is derived from it (Tj is non- vanishing. The term 
X is a parameter which remains arbitrary and that at a special value (x = Xc — 1/6), as we 
have mentioned above, guarantees the renormalizability of the model at leading order in the 
expansion in k. The value x — Xc is usually termed " conformal coupling" , in close analogy to the 
value necessary for the Lagrangian of a gravity-coupled scalar to exhibit conformal symmetry. 
Notice that the Standard Model Lagrangian is not scale invariant, due to the quadratic term of 
the Higgs field in the scalar potential, but would be such if this term were omitted. We will work 
in the almost flat space time limit, in which deviations from the flat metric r]^ v = (+,—,—,—) 
are parametrized in terms of the gravitational coupling k, with k 2 = IQuGn and with Gn being 
the gravitational Newton's constant. At this order the metric is given as g^ u = r]^ u + nh^ u , with 
h^ v denoting the fluctuations of the external graviton. 

The interaction between the gravitational field and matter, at this order, is mediated by 
diagrams containing a single power of the energy momentum tensor (EMT) T^ v and multiple 
fields of the Standard Model. The tree level coupling is summarized by the action 



S,^--J^V, (2) 

where denotes the symmetric and covariantly conserved EMT of the Standard Model La- 
grangian, embedded in a curved space-time background and defined as 

2 5 (S S m + 5 1 / 



(3) 

9=V 
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The complete EMT of the Standard Model, including ghost and gauge-fixing contributions can 
be found in [TT] . 

The fermionic part of the EMT is obtained using the vielbein formalism. Indeed the fermions 
are coupled to gravity by using the spin connection Q induced by the curved metric g^ u . This 
allows the introduction of a derivative which is covariant both under gauge and diffeomor- 
phism transformations. Generically, the Lagrangian for a fermion (/) takes the form: 

C f = % -^{V^) - Z -(V^)^ -mW, (4) 

where the covariant derivative is defined as = <9 M + + with denoting the gauge 
field. The spin connection takes the form 

^ = \° ab v:v bl ,, (5) 

where V is the vielbein, the semicolon denotes the gravitationally covariant derivative and a ab 
are the generators of the Lorentz group in the spinorial representation. The latin indices are 
Lorentz indices of a local free-falling frame, as in Cartan's formulation. 

In establishing the Feynman rules for the perturbative expansion around an almost flat 
background, it is convenient to work in a gauge in which the bilinear mixing terms between 
gauge bosons and their longitudinal Goldstone fields, in the broken EW phase, are removed. 
This induces some modifications respect to the gauge, usually chosen in the computations 
of EW corrections in the flat space-time. In fact, the gauge-fixing Lagrangian in a curved 
gravitational background acquires a new contribution not present in the case of flat space. 
This is due to the promotion of the ordinary flat derivative to a covariant one [17] . which 
requires the addition of the Christoffel connection r° 

d^A v -> X? M = d^A v - T« u A a . (6) 

This is the only term which needs to be varied in the gauge field Lagrangian - together with 
the determinant of the metric y/^g - respect to the background field g^ u . The field strengths, 
for instance, are not affected by the Christoffel connection because of their antisymmetry on 
the Lorentz indices. The gauge-fixing Lagrangian is then given by 




where the gauge-fixing functions are defined as 

p = g ^(d,Wl-T« u W l a ) + l ^{rfa l {H)-{rf)a l H) . = 1,2,3. (8) 
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In the previous equations g, g', and are the coupling constants and the fields of the 
SU{2) W and U(l)y gauge groups respectively, while a 1 are the Pauli matrices and (H) is the 
vacuum expectation value of the Higgs doublet. 

As we have already discussed, the non-minimal coupling of the Higgs scalar, accounted for 
by the Lagrangian Si, generates an extra contributions to the EMT. For this purpose we recall 
that in the broken EW phase, the ordinary parametrizations of the Higgs field 

and of its conjugate W, are expressed in terms of h, and <fi ± , which denote the physical 
Higgs and the Goldstone bosons of the Z and W s respectively, v denotes, as usual, the Higgs 
vacuum expectation value. This expansion generates a non-vanishing EMT, induced by Si, 
given by 

Tj U = -2x(d"d v - rj^n)H^H = -2 X (d>*d u - ^ U U) (y + y + + <T +vhj . (10) 

As we have already mentioned, the most important aspect of the \ — case ^ s the renor- 
malizability of Green's functions with an insertion of EMT and scalar fields on the external 
lines. These are found to be ultraviolet finite only if Tj V is included [211 1221 Ej. Our results, 
however, are presented for an arbitrary \- 



3 The perturbative expansion 

We will be dealing with the Tff (diagonal) fermion case. We introduce the following notation 

= t(p 2 \T^(0)\ Pl ) (11) 

to denote the general structure of the transition amplitude where the initial and final fermion 
states are defined with momenta p± and p 2 respectively. The external fermions are taken on 
mass shell and of equal mass p\ = p\ = m? . We will be using the two linear combinations of 
momenta p = p\ + p 2 and q = p\ — p 2 throughout the paper in order to simplify the structure 
of the final result. 

The tree-level Feynman rules needed for the computation of the T^ u vertex are listed in 
Appendix EH and its expression at Born level is given by 

TT = \u{p 2 ){^p v + Yp»}u( Pl ) . (12) 
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Our analysis will be performed at leading order in the weak coupling expansion, and we will 
define a suitable set of independent tensor amplitudes (and corresponding form factors) to 
parameterize the result. 

The external fermions can be leptons or quarks. In the latter case, since the EMT does not 
carry any non-abelian charge, the color matrix is diagonal and for notational simplicity, will 
not be included. 

We decomposed the full matrix element into six different contribution characterized by the 
SM sectors running in the loop diagrams 

fv> = + + fr> + f •£> + f£ ( 13 ) 

where the subscripts stand respectively for the gluon, the photon, the Higgs, the Z and the W 
bosons and the counterterm contribution. Concerning the last term we postpone a complete 
discussion of the vertex renormalization to a follow-up section. 

As we have already mentioned, we work in the gauge, where the sector of each massive 
gauge boson is always accompanied by the corresponding unphysical longitudinal part. This 
implies that the diagrammatic expansion of T^ v and is characterized by a set of gauge 
boson running in the loops with duplicates obtained by replacing the massive gauge fields with 
their corresponding Goldstones. 

The decomposition in Eq.( JT3l) fully accounts for the SM one- loop corrections to the flavor di- 
agonal EMT matrix element with two external fermions. 

The various diagrammatic contributions appearing in the perturbative expansion are shown in 
FigJU Two of them are characterized by a typical triangle topology, while the others denote 
terms where the insertion of the EMT and the fermion field occur on the same point. The 
computation of these diagrams is rather involved and has been performed in dimensional regu- 
larization using the on-shell renormalization scheme. We have used the standard reduction of 
tensor integrals to a basis of scalar integrals and we have checked explicitly the Ward identity 
coming from the conservation of the EMT, which are crucial to secure the correctness of the 
computation. 



3.1 Tensor decompositions and form factors 

Now we illustrate in more detail the organization of our results. By using symmetry arguments 
and exploiting some consequences of the Ward identities, we have determined a suitable tensor 
basis on which our results are expanded. For massless vector bosons (gluons and photons), 
and for the Higgs field, because of the parity-conserving nature of their interactions, we have 
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Figure 1: The one- loop Feynman diagrams of the graviton fermion vertex. The dashed lines 
can be gluons, photons, Higgs, Z and W bosons or their unphysical longitudinal parts. The 
internal fermion line can be of the same flavor of the external fermions if a neutral boson is 
exchanged in the loop, otherwise, for charged W corrections, it can have different flavor because 
of the CKM matrix. 

decomposed the matrix elements onto a basis of four tensor structures Oy k with four form 
factors fk as 



TT = i^C 2 (N)"£fk(q 2 )u(P2)OP k u( Pl ), 
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. OL 

'An 



k=i 



fc=i 



with the tensor basis defined as 



^^^rn^f h k ^)u{p 2 )0^ k u{p, 

V k=l 



), 



(14) 
(15) 
(16) 



mrf v , 



(17) 



The form factors for the gluon and the photon contributions are identical, the only difference 
relies in the coupling constant and in the charge of the external fermions. The coefficient C 2 (N) 
is the quadratic Casimir in the iV-dimensional fundamental representation, with C 2 = 4/3 
for quarks and zero for leptons. Q denotes the electromagnetic charge and Gp the Fermi 
constant. Moreover, being the fermion-Higgs coupling proportional to the fermion mass, we 
have factorized m 2 in front of the Higgs form factors. Note that Oy 2 _ A are linearly mass 
suppressed, so that only Oy^ survives in the limit of massless external fermions. 
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Coming to the weak sector of our corrections, because of the chiral nature of the Z and W 
interactions, we have to decompose the matrix elements into a more complicated tensor basis 
of six elements as 

*T = i £ fltf) 0^u( Pl ) , (18) 

fc=i 

6 



1 W ~ L 



167T 



7 -^J2^^ u ^°cMPi), (19) 
fc= i 



where we have defined 



= (t^ + tV)^, 

OgS = W + 7> M ) P« , 

= mtT, 

= m (pM g " + ^ p ") 7 5 . (20) 

The most general rank-2 tensor basis that can be built with a metric tensor, two momenta (p 
and q) and matrices 7 M ,7 5 has been given in [11]. The basis given in Eq. (120|) . compared to 
the flavor- changing case, is more compact. We have imposed the symmetry constraints on the 
external fermion states (of equal mass and flavor) and the conservation of the EMT, discussed 
in section HI For the form factors appearing in we introduce the notation 

/iW) = £^/Vf*W.*/) ( 21 ) 
/ 

where Vif is the CKM mixing matrix, with the indices i and / corresponding to the flavor of 
the external and internal-loop fermions respectively, and Xf — mj/m^, where nif and m\y 
stand for the masses of the fermion / and W respectively. 

We have extracted a single mass suppression factor coming from the contribution of the 0^_ & 
operators. In the T^ v matrix element, the leading terms in the small external fermion mass 
are given by the first two form factors. The situation is different for the W case, in which only 
the first form factor is the leading term, being ji suppressed as m 2 . We have decided not to 
factorize the m 2 term in order to make the notations uniform with the Z case. 
We remark that the expressions of the form factors is exact, having kept in the result the 
complete dependence from all the kinematic invariants and from the external and internal 
masses. 
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4 The Ward identity from the conservation of the EMT 



In this section we simply quote the consequences of the conservation of the energy-momentum 
tensor that are contained in the Ward identities satisfied by the matrix elements defined above. 
As widely explained in [17], we can derive a master equation for the effective action T, the 
generating functional of all the 1-particle irreducible (1PI) graphs. We give more details of the 
derivation in Appendix [B] The Ward identities for the various correlators are then obtained via 
functional differentiation. We consider a generic theory with a scalar 0, a gauge field A a and 
a fermion ip. As explained in the Appendix [Bj C can be taken to represents all the scalar and 
ghosts contributions, while A a is a short-hand notation to indicate all the contributions from 
the gauge bosons. 

Imposing the invariance of the generating functional under a diffeomorphism transformation of 
T we have 

2 1 b(p c 5(f)l 5A ca \oA ca 

- s ^ - wf^ + & {§?"*• - ) ■ (22) 

where <r a/3 = [7 Q ,7^]/4 and the subscript c identifies the classical fields. This equation can be 
straightforwardly generalized to the entire spectrum of the SM. 

By a functional differentiation of Eq. (|22p with respect to the fermion fields and after a Fourier 
transform to momentum space we obtain 

qpfr = ^{rfr^-tfr^ (23) 

where Tjj(p) is the fermion two-point function, diagonal in flavor space, given explicitly in 
Appendix [Ql The perturbative test of this equation is of great importance for testing the 
correctness of our results. 

At this point we would like to stress that, as a strong test of our results, we have explicitly 
computed all the form factors entering in the matrix elements of T^ h w z and checked that 
they satisfy the Ward identity in Eq. (T23|) . However, as a consequence of Eq. (|23|) . not all the 
form factors are independent quantities. Therefore, for practical text purposes, we will present 
only the analytical results for the relevant independent subset of form factors, thus reducing 
the number of contributions to the Tff vertex. The other form factors can be derived quite 
straightforwardly by using the Ward identity given above, which determines a set of relations 
among them that will be presented in the next sections. 
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5 Renormalization 



The counterterms needed for the renormalization of the vertex can be obtained by promoting 
the counterterm Lagrangian to the curved background. The counterterm Feynman rules for 
the matrix element with the insertion of the EMT are easily extracted in the usual way and in 
our case, for a chiral fermion, we have 



rp^LV 



*(P2|7ct(0)|;Pi) = -u(p 2 



5z L or 



CI 



SZrO^ + A—O^ 
02 m 



C3 



(24) 



with 5Zl, 5Zr and 5m being the fermion wave function and the mass renormalization constants 
respectively, while 0^[_ 3 are defined in Eq.( l20|) . 

For vector-type interactions, in the gluon, photon and Higgs sector, 5Zl = 5Zr and the expan- 
sion of the T(^p matrix element naturally collapses to only two operators, O v \ = Oq[ + O 



C2 



and 0(^2 



O&of Eq.(H 



We have checked that the renormalization of the parameters of the SM Lagrangian is indeed 
sufficient to cancel all the singularities of the matrix element, as expected. As it can be 
easily seen from Eq. (124j) . the form factors involved in the subtraction of infinities are just the 
first two for the gluon, the photon and the Higgs, and the first three for the massive gauge 
bosons. This is in agreement with simple power counting arguments. 

We have used the on-shell scheme, where the renormalization conditions are fixed in terms of 
the physical parameters of the theory to all orders in the perturbative expansion in the EW 
coupling constants. These are the masses of physical particles, the electric charge and the CKM 
mixing matrix. The renormalization conditions on the fields, which allow the extraction of the 
wave function renormalization constants, are satisfied by requiring a unitary residue of the full 
2-point functions on the physical particle poles. For the fermion renormalization constants we 
obtain the following explicit expressions 



5Z L 
5Z R 
5m 



-ReS L ( 



m 



-ReS 



R (m 2 ) 



2 d ~ 

2 d ~ 
op z 



m~ 
—Re 
2 



E L (m 2 ) 



+ E 



R 'm 2 ) 



+ 2E^ 



) + £V) + 2£V) 

'■) + eV) + 2£V; 
> 2 ) 



(25) 
(26) 
(27) 



where the E L ' R ' S functions are the fermion self-energies defined in Appendix [Dl The symbol Re 
gives the real part of the scalar integrals appearing in the self-energies but it has no effect on 
the CKM matrix elements. If the mixing matrix is real Re can be replaced with Re. 
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6 Form factors for the T^ v matrix element 

6.1 The massless gauge boson contribution 

We give the four form factors for the massless gauge boson cases, namely the gluon and the 
photon contributions. They depend on the kinematic invariant q 2 , the square of the momenta of 
the graviton line, and from the dimensionless ratio y = m 2 /q 2 . The form factors are expressed 
as a combination of one-, two- and three-point scalar integrals, which have been defined in 
Appendix. [Ej and are given by 

2 4j/(2y + l) (8y(7y-4) + 3) , 2) 2(y - 1) 2 , 

h{q) ~ "3(l-4y)2 + WO^WV Ao[m) - 3(1-4^°^ '°' ) 

+ ( 48y"-12 ) g ° ^ ^ + ^ 2(2y " 1)Co (°' m2 ' m ') + 2q [l-^ V C * ^ °' °) ' 
2 4(y(32y - 23) + 3) 2 (32y 2 -26y + 3) 2 2(10?/ -1) , 2 , 

/2(9) = 3(1 -4y) 2 + 3^(1-4^ A ( ffl ) + 3(rW Bo ( g ' Q ' Q ) 

+|^o (g 2 , m 2 , m 2 ) + (1 8 ^ )2 Co (m 2 , 0, 0) , 

2 4y(8y + 19) - 6 4(8y-17) , ( ^ 2(26y + 1) 2 - 

/3(9) = 3 g 2 (4y - 1)3 " 3g*(4 y - 1)3 ^ < m ) ~ 3 g 2 (4y - l)3 g ° ^ » '°) 

+ 3?(%^y o(g ,m ' m )_ (4^iF Co(m '°' 0) ' 

_ 80j/ 2 + 68y _ 9 (20^-1) +3) / ax , (2-20y) , , , 

- ^B (q 2 , m 2 ,m 2 ) - ^^C (m 2 , 0, 0) . (28) 

It is interesting to observe that not all the four form factors are independent because the Ward 
identity imposes relations among them. In fact, specializing Eq. (!23|) to the massless gauge 
bosons contributions, we obtain 

/ 2 (g 2 ) + q 2 h(q 2 ) = ~\ [S 3 L /7 (m 2 ) + Sf /7 (m 2 ) + 2S^(m 2 )] , (29) 

as it can be checked from the explicit expressions given above. This relation can be used to test 
the correctness of our results and to reduce the number of independent form factors. We recall 
once more that the E's denote the fermion self-energies which have been collected in Appendix 

6.2 The Higgs boson contribution 

In this section we present the results for the contribution of a virtual Higgs. As in the previous 
case, they are expanded in terms of scalar integrals and of the kinematic variables q 2 , y = m 2 /q 2 
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and Xh = m 2 /ml, where rrih is Higgs mass. 

As we have already mentioned, the conservation of the EMT induces a Ward identity on the 
correlation functions. This implies a relation between the form factors, which in the Higgs case 
becomes 



££(m 2 ) + ££K) + 2E£(m 2 ) 



^R/^2 



(30) 



Obviously, this equation has the same structure of the Ward identity found in the massless 
gauge bosons case, having expanded 2^ on the same tensor basis of Eq. (TT7j) . 
Notice that the f£ and f% form factors depend on the x parameter of Sj. This is expected, 
because the Higgs field can also couple to gravity with the EMT of improvement Tj V defined 
in Eq. fllOp . The Feynman rules for a graviton-two Higgs vertex are then modified with the 
inclusion of the \ dependence and affect the diagram represented in Figfjjb), where this vertex 
appears. We obtain 



h(„2\ 



3x h -8y- Ax h y 



+ 



12x h (l-4y) 3q 2 (l-Ay) 



A (m 2 h ) - Ao (m 2 ) 



+ 



+ 



+ 



12x2(l_ 4 y)2 
1 

6x2(1 -4y)2 



x\ + 8(x h (2Qx h + 3) - 3)y 2 - 2(28x h + 3)x h y 



K2 I 2 2 2\ 

B [q ,m ,m 



x\ + 4(3 - 8x h )y 2 + 4(2x h - l)x h y 



B (q ,m h ,m h ) 



2x h (\-Ay) 2 



x h (A - 32y) + 12y + 1 



B (m ,m ,m h ) 



q 2 y 



+ A(x h (Ax h - 3)(4x h + 1) + l)y 2 + (8(1 - 4x h )x h + 3)x h y 
q 2 y{xh - 2y) 



2x3(1 -4y) 2 
C (m 2 h ,m 2 ,m 2 ) 



Axi 



+ 



h/„2\ 



x\(l - Ay) 2 
40x^y - 9x h - Ay 



x 2 ~ ^x h y + y) Co (m 2 ,m 2 h , m 2 h ) , 
4 



+ 



3x h (l- Ay) 3q 2 (l-Ay) 

Bo (q 2 ,m 2 ,m 2 ) + 
2 



+ 9x h y(l - Ay) + 6y 2 
x B (q 2 ,m 2 h ,m 2 h ) + 



Ao (ml) - Ao (m 2 ) 
1 



3x2(1 -Ayf 



2x 2 h (l -Ayf 



x h (l - Ay) 2 
Aq 2 y(-Ax h y + x h + y) 2 



3x2 (1 _ 4y)2 
4(1 - Ax h )y 2 + Q Xh y -x h + y 



x z h + 4(3 - 2Qx h )y* + 8(x h + l)x h y 



Bq (m 2 ,m 2 ,ml) 



x\(l - Ay) 2 
q 2 (x h (8y - 1) - 2y) 



x'Hl-Ay) 2 



C (m h ,m ,m ) 



x 2 h (2y - 1) + 8x h y 2 - 2y 2 



+ X 



I -Ay 



Bo{m 2 ,m 2 ,m 2 h )-Bo{q\mlml) 
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3q 2 x h (l - Ayf 



+ 



3qHl - Ay) 2 y 
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+ 
+ 
+ 
+ 
+ 
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3q 2 x h (l-Ay) 3g 4 y(l - 4y) 
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n ( 2 2 2\ 
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+ 



B (q 2 ,m 2 ,m 2 ) 



1 



+ 8(x h + l)x h y 



3q 2 x 2 h {l - Ay) 2 
1 



3q 2 x 2 h {\ - Ay) 2 
x\ + 4(3 - 2Kx h )y 2 



24(1"%)' 



q 2 x h (l - Ay) 2 



V U (m 2 m 2 rr> 2 \ J- + V ^ X hU) 2 p 1 2 2 2\ , 

x £> [m ,m ,m h ) H 3— — — ^ C [m h ,m ,m ) + 



x h (A{5 - 8y)y - 2) + 2y(4y - 5) + 1 
(x h (8y - 1) - 2y) 
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C (m 2 ,m 2 h ,m 2 h ) +\ 
A(x h + 2y- 8x h y) 



x 3 h {l-Ay) 2 



4(2y-i) 



9 2 (1 -%) 



£ (q , m h , m h ) - B (m ,m , m h ) 



•n I 2 2 2^ 

-C (m ,m h ,m h 



x h (l-Ay) 

6.3 The Z gauge boson contribution 



(31) 



Coming to the form factors for the Z boson contribution, which are part of T% v ', these are 
given in terms of the variables q 2 , y = m 2 /q 2 and Xz = m 2 /m|, with the parameters v and a 
denoting the vector and axial-vector Z-fermion couplings. In particular we have 



v = I 3 - 2s 2 w Q , a = h 



2 2,2 

c = v + a 



(32) 



where ^3 and Q are, respectively, the third component of isospin and the electric charge of the 
external fermions, while sw is the sine of the weak angle. 

In this case, the structure of the Ward identity is more involved than the previous cases, 
being T^ v expanded on a more complicated tensor basis, Eq.( l20i) . We obtain two relations 
among the form factors that we have tested on our explicit computation, which are given by 



/a 
h 



/f + ^ + ^Sf(m 2 )-S|(m 2 ) 



E|(m 2 ) + E| (m 2 ) + 2£|(m 2 
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(33) 
(34) 



Also in this case we have a dependence of the result on the parameter x, which appears in f§ 
and hence in As for the Higgs field, also the gravitational coupling of the Z Goldstone 
boson acquires a new contribution coming from the term of improvement Tj, shown by the 
Feynman diagram in Fig{T](b). 

Here we present a list of the explicit expressions of /-f , ff, f§ and f$ while f£ and f£ can 
be obtained using the Ward identity constraints of Eq. (l33|) and Eq. flMl) . We obtain 
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+ 



+ 



q 2 y 
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Aq (m 2 z ) - Ao (m 2 ) 
x z {{Ay - l)x z {-Ay (8a 2 + Aav + llv 2 ) + 2a 2 (4y -l)x z + 17 (a + v ) 2 ) 
+ 6y (4j/ (5a 2 + 8au + 7u 2 ) - 7(a + u ) 2 )) - 24y 2 (a + v) 2 B (g 2 , m 2 , m 2 ) 

2a 2 i/ r 

x z {x z {-2y (15a 2 + 20av + v 2 ) + a 2 (8y + l)x z + 64a 2 y 2 + 2(a + v) 2 ) 

So (g 2 ,m|,m|) 
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q 2 y 



(1 - Ay) 2 x 2 z 
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S (m ,m ,m z ) 



Aq 4 y 2 



(1 - Ay) 2 x% 
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/i / 2 2 2 \ 
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91/ 
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Aq (m 2 ) - A {m 2 z ) 
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Mj 
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+ 2a 2 (64y 2 - 3Ay - 3) + c 2 (26y + 1))) - 30c 2 y 2 B (q 2 ,m 2 z , m|) 



+ 



Ay 
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2y 
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((8y - 9)x z - 8y) + 
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3(4y - 1)4 

2avy 
3(\-Ay) 2 x\ 
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3q 2 (\-Ay)x z 
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2avy 



+ 



(1 - Ayfx 2 z 
16avq 2 y 3 
(l-4y) 2 4 
Aav q 2 y 



B (m 2 ,m 2 ,m 2 z ) 



(8y + 2)x z -2y + l 
x z (x z {Ax z - 3) - Ay + 1) + y 



C (m 2 ,m 2 z ,m 2 z ) 



Ay) 2 xt 



{{Ay - l)x z - y) 



{8y - 5)yx z + {Ay - \)x\ - Ay 2 



Co{m 2 z , m 2 , m 2 ) .(35) 



6.4 The W gauge boson contribution 



Finally we collect here the results for the matrix element. They are expressed in terms of 
scalar integrals and of the kinematic invariants q 2 , y = m 2 /q 2 , xw = m 2 /m 2 v and Xf = m 2 /m 2 , 
where mj is the mass of the fermion of flavor / running in the loop. 

As in the Z boson case the conservation equation for the EMT implies the following relations 
among the form factors 



w 

2 

w 



fY + ffe W + \ 



2 f W 
1 



J W 



m 



m 



+ 2S^(m 2 ) 



(36) 
(37) 



which we have tested explicitly. Also in this for the form factors with the exchange of 

a Z boson, , and hence f^, depends on the parameter x- 

We recall that the 0^3_ 6 operators are characterized by a linear mass suppression in the 
limit of small external fermion masses, while Oq^i even if not explicitly shown, has a quadratic 
suppression, which is present only in the W case. Therefore, the leading contribution, in the 
limit of external massless fermions, is then given by the first form factor fY alone. 



We present the explicit results for the and to F^ in Appendix |A], while F^ and 
can be computed using the Ward identities of Eq. fl36|) and fl37|) . The form factors f^ 
are obtained from F^ multiplying by the CKM matrix elements and then summing over the 
fermion flavors, as explained in Eq.( l2TT) . 



F w 



7 Infrared singularities and soft bremsstrahlung 

Here we provide a simple prove of the infrared safety of the Tff vertex against soft radiative 
corrections and emissions of massless gauge bosons. 

An infrared divergence comes from the topology diagram depicted in Figfjja) with a virtual 
massless gauge boson exchanged between the two fermion lines and it is contained in the three- 
point scalar integral Cq(0, m 2 , m 2 ). If we regularize the infrared singularity with a small photon 
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(or gluon) mass A the divergent part of the scalar integral becomes 

A 



C o (0,m 2 ,m 2 



J- s 



■ft 



{ - 2 log — log x s + . . . ) 



(38) 



where the dots stand for the finite terms not proportional to log^, and x s = —j^ with 
= a/1 - 4m 2 /q 2 . 

In the photon case the infrared singular part of the matrix element is then given by 



rpfiU 

7 



47T 



a 

47T 



y 



<2y-- 

-(2J/ -1 



log — log x s 

m 



x s A 
log — log x t 

1 — x z s m 



u(p 2 ) Oy>(pi) + . . . 

1 "I"---, 



(39) 



which is manifestly proportional to the tree level vertex. On the other hand, the gluon contri- 
bution is easily obtained from the previous equation by replacing aQ 2 with a s C2(N). 

For the massless gauge boson contributions there is another infrared divergence coming from 
the renormalization counterterm. Its origin is in the field renormalization constants of charged 
particles arising from photonic or gluonic corrections to the fermion self energies. For example, 
in the photon case we have 



5Z r 



IR 



5Z 



R 



IR 



Air L m 



+ ... 



(40) 



However the processes described by the T^ v matrix element alone are not of direct physical 
relevance, since they cannot be distinguished experimentally from those involving the emission 
of soft massless gauge bosons. Adding incoherently the cross sections of all the different pro- 
cesses with arbitrary numbers of emitted soft photons (or gluons) all the infrared divergences 
are expected to cancel, as in an ordinary gauge theory [23l EH [25] . This cancellation takes 
place between the virtual and the real bremsstrahlung corrections, and is valid to all orders in 
perturbation theory. In our case one has to consider only radiation of a single massless gauge 
boson with energy ko < AE, smaller than a given cutoff parameter. 

For definiteness we consider the emission of a photon from the two external fermion legs. 
The gluon already mentioned, is easily obtained from the final result with the replace- 

ment aQ 2 — > a s C2{N). In the soft photon approximation the real emission matrix element, 
corresponding to the sum of the two diagrams depicted in Fig. [2] is given by 

-e(k)-pi e(k)-p 2 ~ 



M S oft = M (eQ) 



(41) 



k-pi k-p 2 

where k and e(k) are the photon momentum and polarization vector respectively, while the sign 
difference between the two eikonal factors in the square brackets is due to the different fermion 
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(a) (b) 
Figure 2: Real emission diagrams of a massless gauge boson with momentum k. 



charge flow of the diagrams. Here A4q is the Born amplitude which factorizes, in our case, as 

M = A, u fr, (42) 

where Tq V is the tree level graviton vertex defined in Eq. (fl2|) and A^ v is the remaining amplitude 
which does not participate to the soft photon emissions. 

The cancellation of the infrared singularities occurs at the cross section level, therefore we 
have to square the soft photon matrix element, sum over the photon polarization and integrate 
over the soft photon phase space 



do 



soft 



a o / d 2 k 
2n z 



|fc|<A£ 



Pi 



+ 



Pi 



- 2- 



Pl -P2 



(k ■ pi) 2 (k ■ pi) 2 k ■ p\k ■ p 2 - 



(43) 



where the infrared divergence is regularized by the photon mass A which appears through the 
photon energy k = \J\k\ + A 2 . 
The generic soft integral 



Hi 



d 2 k 2pi ■ pj 
\< AE 2k k-p i k-p j 



(44) 



has been worked out explicitly in |26j, here we give only the infrared divergent parts needed in 
our case 



hi 

h2 



AE 

7 22 = 47T log — h 



A 



l7T 1 



3-> R 



2m Vl-^ 0g T l0gXs + - 



(45) 



with x s = — and f3 = a/1 — Am 2 /q 2 . 

Using the previous results we obtain the infrared singular part of the soft cross section 



, a ^2 f , A£ 8 . . x s , A£/' , i 

dcr sof t = -da —Q <8\og — — + -(2y - 1)- log — — \ogx s + ... \ 

47r I Ay 1 — x{ A J 



AE 



(46) 
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where dots stand for finite terms. 

Exploiting the fact that the infrared divergences in the one-loop corrections and in the 
counterterm diagram multiply the tree level graviton vertex Tq V , so that they contribute only 
with a term proportional to the Born cross section, we obtain 



d(T V i r t + docT + da 



T soft 



dcro j-Q 2 

4:71 



1 + -(22/- 1) 

y 



X a 



1 — x\ 



■ log x s 8 log 



m 



+ ... 



(47) 



for the photon case and an analog result for the gluon contribution. The sum of the renormalized 
virtual corrections with the real emission contributions is then finite in the limit A — > 0. 



8 Conclusions 

We have computed the one-loop EW and strong corrections to the flavor diagonal graviton- 
fermion vertices in the Standard Model. This work is an extension, to the flavor diagonal case, 
of previous related study in which only the flavor-changing fermion graviton interactions had 
been investigated. The result of our computation has been expressed in terms of a certain 
numbers of on-shell form factors, which have been given at leading order in the EW expansion 
and by retaining the exact dependence on the fermion masses. We have also included in our 
analysis the contribution of a non-minimally coupled Higgs sector, with an arbitrary value 
of the coupling parameter. All these results can be easily extended to theories with fermion 
couplings to massive graviton, graviscalar and dilaton fields. 

Moreover, we proved the infrared safety of the fermion-graviton vertices against radiative 
corrections of soft photons and gluons, where the ordinary cancellation mechanism between 
the virtual and real bremsstrahlung corrections have been generalized to the fermion-graviton 
interactions. 

There are several phenomenological implications of this study that one could consider. 
Beside the possible applications to models with a low gravity scale, which would make the 
corrections discussed here far more sizeable, one could consider, for instance, the specialization 
of our results to the neutrino sector, a definitely appealing argument on the cosmological side. 
Another possible extension would be to include, as a gravitational background, also a dilaton 
field, generated, for instance, from metric compactifications. We plan to return to some of these 
open issues in the future. 
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A Appendix. The structure of the W-boson form factors 

We collect in this appendix the expression of the form factors generated by the exchange of a 
W-boson in the loop. They are given by 
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+ 4yx w 



+ 



(8y - 3) (x f (x w + 2) + x w ) 



A ( m w) ~ A) (w*/) 



3g 2 (l - 4y)x/x w 

3 (x^ + 2) (l2y 2 (x w - l) 2 - 4y (xw - 3) x w + x 2 ^ - x 2 x w 



3(1 - 4y) 2 x 3 x 3 w 

x (8(2j/ + l)yx w + (4y(lly - 7) + 5)x 2 w + 36y 2 ) + 4(2 - lly)yx f x 3 w + 12y 2 x^ 



2y 



x ^( g , m/>ro/ ) + -^— 



^ 0*7 0*7 ((y(% + 5) - 2)x/ + 2(10 - 9y)y 



- 5) + 9y(2y - 1)) - 6y 2 ) + (4y - l)x^x w ((7?/ - 8)x/ + y) + 2yx 2 x w ((13?/ 



l)x/ + 9y)-12y 2 x^ 



x 2 , (% ((1 - 2y(4y 



(1 - 4y) 2 x 2 x 2 w 

+ l))x w + 2(9 - 4y)y - 5) + 4y(4y - 5) + 2) + x/x w (4(1 - 2yfx w + 2(5 - 4y)y 

2q 2 y 



- l) + (2(5-4y)y-l)x w 



S (m 2 ,m^m^)- 



(1 — 4y) 2 x 4 x 4 



10y 3 (x/ + 



- 2yxfx^ ((2y + l)yx f + (y(3y + 4) - l)x} + 3y 2 ) - y (x f + 1) x f x 3 w (2(1 

- 2y)yx / + (2(y - 6)y + 3)x 2 + 2y 2 ) + x A w (2(1 - 2y)yx/ + (2(y - 2)y + l)x) 

+ 2y 2 ) (x f (yx f - 2y + 1) + y) - 4y 3 xj C (m 2 ,m 2 w , m 2 w ) - ] q -| 4 4 : 

j v 1 4 yj x f x w ~ 



2xfXw 
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x ((y - l) Xf + y) + (x f - 1) x 2 w {{2y - l)x f - 2y) - Ayx) 



- 2yx f x w (xw + 1) + yx"w 



/ (y(xw - i) 2 + x w ^j 

C (m w , mj, m 2 f ) + x| ^ ^4^)^ B (q 2 ,m 2 w ,m 2 w ) 



Bq (nr. 1 



+ 



8q y 



{Ay - l)x 2 f x w '_ 



yx f (xf + 1) - x w (xf (yxf -2y + \) + y) 



x C (m},m 2 w ,m 2 w ) 



W (2\ 



y 



6(1 — Ay)x\x\ 



+ 



(8y - 3) (xf (x w ~ 2) - x w ) 



2x f x w ((Ay -9)x f + Ay) + (x f -l)x 2 w (3x f + 8y) + 16yx 2 f 



6q 2 (l — Ay)xfXw 



Aq (m w ) - Ao (m 2 ) 



6(1 - Ay^xjx^ 



(x f - 1) 



X X 



3 W (2(5 - 8y)yx f + (2y(12y - 7) - l)x) + 2Ay 2 ) + 2x 2 f x 2 w (((53 - A8y)y - 17)x f 



+ y(16y + 23)) + 12yxjx w ((Wy -7)x f + 6y) - A8y 2 xj 



Bo (g 2 , 



2 _2 



m f ,m f 



) 



3(1 - Ay^xjx^ [ 



(xf-1) x 3 w (12(1 - 2y)yx f + (Ay(3y - 1) + l)x) + 12y 2 ) 



+ 8x 2 x^ ((4(i/ - 2)y + l)x/ - 2y(2y + 1)) - 4yx^x w (7(y - l)x/ + 9y) + 24y 2 x) 



i3 (^TO^m 2 ^) 



2(1 - 4y)2x 2 x 2 y 
+ z/^w (~4yx w + 2y - 1) + (2y - ljx 2 ^ 



x 2 (x w (2y (x w + 7) + x w + 5) - Ay + 2) 

2?V 



£0 ( 



m , m 



,m w ) - 



(1 — Ay) 2 x f x- 



w 



2yx} 



x x w ((3y - 2)x f + 5y) + yx 2 x^ (x/ ((2y + 7)x f + 4y + 5) - 6y) - 2yx f x w (x f (x f 
x ((3y + l)x f - 5y + 5) + y + 2) + y) + (x f - I) x$y ((2y - l)x/ - 2y) (x/ (yx/ 



2y + 1) + y) - Ay l x) 



C (m 2 ,77^,771^) 



q 2 y 



2y 2 xjx w ((Uy-9)x f 



(1 - Ay) 2 xjx^ 

+ Wy) - yx 2 f x 2 w ((16y - 25)yx f + 3(y(12y - 13) + 4)x^ + 12y 2 ) - y (x f - l) 2 x A w 
x (-yx f + (y(4y - 3) + l)x) - Ay 2 ) + 2x/x^ (x f (x f ((y(2y(5y - 6) + 5) - l)x/ 



- 6y 3 + 3y) -2y 2 (y + 2)) - 2t/ 3 ) - 8y 3 x) 



,m f ,m f 



)• 



(48) 



B Appendix. The Ward identity 

In this appendix we fill-in some of the gaps in the derivation of the Ward identity satisfied by 
the effective action of the Standard Model. In order to simplify our notations, as mentioned in 
the main section, we take as an illustration a theory with a spin-1 gauge field (A^), a charged 
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scalar (0) and a fermion (■?/>). In the case of the Standard Model the extension of this analysis is 
straightforward but quite lengthy. We have replicas of the spin-l's (i.e. A® = (A g , A^,W\ Z)) 
on which we implicitly sum over, while the scalar summarizes both the Higgs and the ghost 
contributions (0 = ("H,u/), (<jr = (H^ , u) 1 )) and ip all the fermions. 
In these condensed notations, the conservation equation of the EMT of the Standard Model 
takes the following off-shell form 

SS n , -5S 1 _„ f 8S , - 8S\ .„ 8S 



"" 8^ ur 5*p ' 2 V# M <W ^ 



+ s "l*S)-ffw-w§ < 49 > 



where we are implicitly summing on all the spin-1, scalars and fermions. The off-shell relation 
( H9|) can be inserted into the functional integral in order to derive some of the Ward identities 
satisfied by the Tff correlator. Whence we define the generating functional of the theory in 
the fluctuations of the background gravitational metric (h^ u ) 

Z[J,r],rjXX\h^] = J VAT>ipVipV(j)Vtf exp {i J d A x (C SM + 

+^ + ^ + C t0 + 0t c + } . (50) 

The EMT, obviously, is chosen to be the symmetric one and on-shell conserved. We have 
denoted with J,r],fjX,( the sources of the gauge field(s), the fermion and antifermion fields 
and scalar and its conjugate respectively. The generating functional W of the connected Green's 
functions is, as usual, denoted by 

exp i W[J, V, V, C, C\ = — ( 51 ) 

(normalized to the vacuum functional). The effective action, defined as the generating func- 
tional T of the 1-particle irreducible amplitudes is obtained from W by a Legendre transfor- 
mation with respect to all the sources, except, in our case, the metric fluctuation h^, which is 
taken as a background external field 

T[A C , 4, Vc € (Pc V] = W[J, v , v, C, C f , V] " / d ' x ( J » A c + #o + + CVc + <f>tC) • (52) 



As usual, we eliminate the source fields from the right hand side of Eq. ( 1521) inverting the 
relations 

K - 83 > ^ ~ 5fj ' ^ ~ 8rj ' 0C " 80 ' 0C " 8( (53) 
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so that the functional derivatives of the effective action T with respect to its independent 
variables are 

5T 5T 5T 6T , 5T 



5A£ <fy> c " <fy c " Scf> c 54I 



and for the source we have instead 



& - w- < 55 > 

The conservation of the EMT given by Eq. (1491) is rewritten in terms of classical fields and 
then re-expressed in functional form by differentiating W with respect to h^. We use Eq. (1491) 
under the functional integral. We obtain 

+ 5W 5W 

+ C^^ + ^C, (56) 

and finally, for the one particle irreducible generating functional, this gives Eq. f[2"2"j) . after using 
Eq. flM} with the help of Eqs. (135]), (}5l|), (155]). 



C Feynman rules 

We collect here all the Feynman rules involving a graviton that have been used in this work. 
All the momenta are incoming 

• graviton - gauge boson - gauge boson vertex 




(57) 

where V stands for the vector gauge bosons g, 7, Z and W. 
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graviton - fermion - fermion vertex 




-t^{r (h + k 2 y + y (h + k 2 y -2rr (h + h- 2m f) } 



(58) 



graviton - scalar - scalar vertex 
5* 



k2- 



k lp k 2a C^-M 2 sV ^ 

i^2 X [(h + k 2 y(h + k 2 y - rr(h + k 2 f] 



(59) 

where S stands for the Higgs H and the Goldstones and ± . The first line is the 
contribution coming from the minimal energy-momentum tensor while the second is due 
to the improvement term. 

graviton - scalar - fermion - fermion vertex 




where the coefficients are defined as 



^hipip ~ 

c L - 

<p Iplp 



e m 
2sw mw 



m. 



V., 



V2s w m w 
e mj 

-i—i= - 

y/2s w m w 



V; 



ipip ' 



c 



e m 

> 

2sw mw 



2h 



c5 



n™ - — 1 
U f~^~ l V2. 



\/2s w m w 

e 



s w m w 



V- 



(60) 



(61) 
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• graviton - gauge boson - fermion - fermion vertex 




7/3 



with 



(JL_ . 

gipip 

r L 



C R - 

gtptp 



ig s T\ 



c 



2S\yCw 

e 



[v + a), C 



C L - 



R 



C R - 



W+ipip 



= I 



V2 



-v., 



s w 



Iplp ? 



2S\yC\y 
i—^r V; 



ieQ , 
v-a), 



y/2. 



sw 



n R - — n 



and v — I 3 - 2s w Q, a = I 3 . 



(62) 



(63) 



The tensor structures C, D and E which appear in the Feynman rules defined above are given 
by 

Cfn/pa Vpp Vvcr Vpcr T]i/p Vpv Vp& ) 

D^paih , k 2 ) = T}^ h a k 2 p - rf G k\k p 2 + ri w k la k 2v - r) pa h M k 2 „ + (// v) 



E^paih, k 2 ) = r]p, v (k lp k la + k 2p k 2a + k lp k 2a ) - r] va k 1/x k lp + r] up k 2 pk 2a + (// -B- v 



(64) 



D Fermion self-energy 

The one-loop fermion two-point function, diagonal in the flavor space, is defined as 
T ff (p) = i [j>P L S V) + i>Pn £ V) + m £ V)] 



(65) 



where the three components S x (p 2 ), with X = L,R, S, are eventually given by the gluon, the 
photon, the Higgs, the Z and the W contributions 



m 2 ^(p 2 ) + ^z(p 2 ) + ^w(p 2 ) ■ (66) 
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The E (p ) coefficients of the fermion self-energies are explicitly given by 



£ 9 V) = Sf (p 2 ) = £^(p 2 ) = £* (p 2 ) = -2S a (p 2 ,m 2 ,0) - 1 
£ 9 V) = £?(p 2 ) = -4 So (p 2 ,m 2 , 0) + 2 , 
££(p 2 ) = X£(p 2 ) = -2Si(p 2 ,m 2 ,m 2 ) , 
£f(p 2 ) = 2S (p 2 ,m 2 ,m 2 ) , 

^(p 2 ) 
^(p 2 ) 



" 4 E V *f V fi [ K + 2m w) Bi (p 2 , m% m 2 w ) + 



ir 



-Am 2 Y,V* f V fi B 1 (p 2 



2 2 \ 

m f ,m w ) , 



£§(p 2 ) 

£|(p 2 ) 

£|(p 2 ) 



-2m%(v + a) 2 



2B 1 (p 2 



,m ,m Z/ 

2 2 2 

- ,m z 



) + 1 - 2m 2 Bi (p 2 



m , m 



2m|(w - a) 2 2 Si (p : 

4 i3 (p~,m~,m z 



-2ml(v 2 - a 2 



) + 1 - 2m 2 Si (p 
)-2 



m , m 



2 _2 _2 

2m 2 So (p~, //'". //' 



(67) 



where v and a are the vector and axial- vector Z-fermion couplings defined in Eq. fl32l) and 



Si (p 2 ,m^,m 2 ) 



m\ — m, 



2p 2 



B (p 2 ,ml,m 2 1 ) - So(0,mo,m 2 ) - ^S (p 2 , mg, m\) . (68) 



E Scalar integrals 

In this Appendix we collect the definitions of the scalar integrals appearing in the computation 
of the matrix element. One-, two- and three- point functions are denoted respectively as Ao, 
S and C with 



#o(p 2 ,™ 2 ,m 2 ) : 



m, 



2 • 




d n l- 



1 



Co(pI, (Pi -P2) 2 ,pl,ml,m 2 1 ,m 2 2 ) = — i d n l 

IK 1 J 



in-./ (I 2 — ml) ((I + pi) 2 — m\) ' 

1 



(I 2 -m^)((/ + pi) 2 -m 2 1 )((l + p 2 ) 2 -ml) ' 

(69) 



Because the kinematic invariants on the external states of our computation are fixed, q 2 = 
{p\ — P2) 2 , Pi = pl = m 2 , we have defined the shorter notation for the three-point scalar 
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integrals 



n I 2 2 2\ n I 2 2 2 2 2 2\ 

C (m , m l5 m 2 ) = C {m ,q ,m ,m , m lt m 2 ) , (70) 
with the first three variables omitted. 
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